Experimental Quantum Imaging exploiting multi-mode spatial correlation of twin 

beams 
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Properties of quantum states have disclosed new and revolutionary technologies, ranging from 
quantum information to quantum imaging. This last field is addressed to overcome limits of classical 
imaging by exploiting specific properties of quantum states of light. One of the most interesting 
proposed scheme exploits spatial quantum correlations between twin beams for realizing sub-shot- 
noise imaging of weak absorbing objects, leading ideally to a noise-free imaging. Here we discuss in 
detail the experimental realization of this scheme, showing its capability to reach a larger signal to 
noise ratio with respect to classical imaging methods and, therefore, its interest for future practical 
applications. 

PACS numbers: 42.50.Ar, 42.50.Dv, 42.50.Lc,03.65.Wj 



I. INTRODUCTION 



The possibility of manipulating quantum states as 
atoms, photons, etc. has recently fostered the develop- 
ing of quantum technologies [ll-[l8| , with very promising 
opportunities for future widespread applications. Among 
quantum technologies the ones based on photons, in par- 
ticularly exploiting entanglement as a resource, have pro- 
duced the most impressive results, such as quantum cryp- 
tography [2|, Q? teleportation [j-l6|, linear-optics quan- 
tum computation 0, entanglement distillation [8[ and 
purification |9l-[lH, quantum metrology [I^- The use of 
quantum properties of light for overcoming limits of clas- 
sical imaging is one of the most promising |12l-[l^. In 
particular it was suggested that, by exploiting the strong 
correlation in noise of entangled modes of light produced 
by Parametric Down Conversion (PDC), the image of a 
weak absorbing object in one branch, eventually previ- 
ously hidden in the noise, can be restored by subtracting 
the spatial noise pattern measured in the other branch, a 
technique that has been christened Sub Shot Noise Quan- 
tum Imaging (SSNQI) [23|. When operating under shot 
noise regime (avoiding the limits of the need of back- 
ground subtraction [2^ or single spatial mode operation 
[l9l-[2ll, l3ll|) the method allows in principle a full recon- 
struction of the absorption pattern of an object with a 
sensitivity superior to that available with classical tech- 
niques at the same illumination level. In order to reach 
this goal in view of important practical applications one 
should reach a very high level of spatial quantum correla- 
tion. In [23| we have demonstrated that this level can be 
effectively reached and, then, in [25| we have presented a 
proof of principle of the sub shot noise quantum imaging 
scheme. Here, in order to provide all the necessary in- 
formation, we present further data, discussing the details 
of the experimental realization. In particular, we show 
how to reach a degree of spatial correlation needed for 
SSNQI performances better than the corresponding clas- 
sical scheme, i.e. a differential classical imaging (DCI) 
scheme in which, instead of the twin beams of PDC, two 
classically correlated beam are used. This achievement 



demonstrates that an experimental setup for SSNQ imag- 
ing is feasible and the practical application of the method 
can be upcoming. In SecHIlwe provide a basic introduc- 
tion to the multi-mode spatial squeezing and the idea of 
SSNQI exploiting PDC, while Sec lIIII describes the exper- 
imental set-up and the spatial sub-shot-noise measure- 
ments. SecHVl considers some effects that spoil the spa- 
tial sub-shot-noise regime and their possible corrections. 
Finally, in SecIVlwe present our quantum imaging exper- 
iment, and two different figures of merit that compare 
the performances of SSNQI and DCI. Some conclusion 
and remarks can be found in Sec JVIl 



II. THE SCHEME 

The scheme 0| for achieving quantum imaging of 
weak absorbing objects, is based on considering two spa- 
tially correlated areas of Parametric Down Conversion 
S'DC) light. This is a non-linear optical phenomenon 
[iSj where a pump beam photon " decays" in two lower 
frequencies photons, usually called signal (s) and idler 
(i), strongly correlated in frequency and momentum since 
they have to conserve the energy and momentum of the 
original pump photons. In particular, the conservation 
of the transverse (to the propagation direction) compo- 
nent of the momentum traduces in a correlation between 
the photon number A/'^(x) and 7Vs(— x) ideally detected 
in any pairs of symmetrical positions x and — x in the far 
field plane (note: we restrict to wavelength near degen- 
eracy wavelength Ag = A^). The degree of correlation is 
quantified by the Noise Reduction Factor (NRF) 
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where (^^O) = (O^) — (O)^ is the mean square fluctuation 
of a generic operator O. The normahzation in Eq. ([1]) 
represents the Shot Noise Limit (SNL) (A^g + iV^), deflned 
as the level of fluctuation associated to the difference be- 
tween coherent beams. In terms of the transmission rj 
of the optical channel (including the quantum efficiency 
of the detector) the theoretical value of NRF for PDC 
is a = 1 — T] (considering for simplicity balanced losses 
Vs = Vi = V that leads to symmetrical statistical proper- 
ties of the detected beams). Therefore, in an ideal case in 
which 77 ^ 1, <j approaches zero. On the other side, for 
the subtraction of two classical beams the degree of cor- 
relation is bounded by adass ^ 1, where the lowest limit 
o'ciass = 1 is reached for coherent beams, or correlated 
beams, generated by a single classical beam separated by 
a balanced beam splitter. In this case the subtraction of 
the detected intensities allows to eliminate the classical 
excess noise contained in the source, while the shot noise 
survives, leading to {5'^{Ni - Ns)) = "^{Ns). The idea 
of sub-shot noise quantum imaging (SSNQI), depicted in 
Fig. [TJ consists in measuring the intensity pattern on one 
branch (e.g. the signal one), where the object has been 
inserted, and then subtracting the correlated noise pat- 
tern measured in the other branch (the idler one) that 
does not interact with the object. The number of photon 
detected in presence of the object in the signal region is 
(A^^(x)) = [1 — a(x)](A^s) where a(x) is the absorption 
in the position x. 

Therefore, in the SSNQI scheme we evaluate the ab- 
sorbtion as 



a(x) = ([iV,(-x)-iV:(x)])/(iV,) 



(2) 



where implicitly we always consider (Ni) = {Ng). 

The Signal to Noise Ratio (SNR) is improved because 
the spatial noise that affects the single beam is washed 
away by the subtraction. In Ref.[lJ the theoretical capa- 
bility of the described method exploiting PDC has been 
compared with the ones of the corresponding differential 
classical scheme in which a coherent beam is split by a 
50% beam splitter. It was shown that this leads to a ra- 
tio between the SNR in quantum and in the differential 
classical imaging of 
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for the same number of photons. We have introduced the 
excess noise En = {{S'^Ni) - {Ni))/{Ni) of the PDC, i.e. 
the noise that exceeds the SNL, also referred to as stan- 
dard quantum limit. For multi-thermal statistics of the 
single branch of PDC light, the excess noise is roughly 
given by {N)/A4, with A4 the number of modes. Equa- 
tion [3] shows that, when the excess noise is negligible (in 
particular a'^E^ <C 1), the SSNQI presents an advan- 
tage respect to a classical differential imaging for a weak 
absorbing object {a -^ 0) as soon as a < 1. This for- 
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FIG. 1: Experimental set up. A 355 nm laser beam pumps 
a type II BBO crystal producing PDC. After eliminating UV 
beam, one correlated beam cross a weak absorbing object and 
it is addressed to a CCD array. The other beam (reference) is 
directly addressed to another area of CCD camera. The total 
transmittance of the optical path is evaluated to 70%. 



mula can be generalized to the case of multi-thermal su- 
per Poissonian source for the differential classical scheme 
(instead of only considering coherent light) just including 
in the numerator under square root the term E^^^^a^ ^ 
that represents the contribution due to the excess noise 
^ciass q£ ^Yie classical light. 

Furthermore, it can be demonstrated that for a < 0.5 
the quantum imaging technique is more advantageous 
than an ideal direct classical scheme, and therefore it 
can be a useful resource for practical applications |26j . 



III. THE EXPERIMENTAL REALIZATION OF 
THE SUB-SHOT-NOISE 

In our experimental setup a type II BBO non-linear 
crystal (1 = 1 mm) is pumped by the third harmonic 
(355 nm) of a Q-switched Nd:Yag laser. The pulses have 
a duration of Tpump = 5 ns with a repetition rate of 10 
Hz and a maximum energy, at the selected wavelength, 
of about 200 mJ. This choice for the pulse duration is 
motivated by the indications of Ref. [1^ for realizing a 
set-up suited for a test of SSNQI, overcoming the limits 
of Ref. |2J]. After eliminating with a prism the resid- 
ual of first and second harmonic components, the pump 
beam crosses a spatial filter (a lens with a focal length 
of 50 cm and a diamond pin-hole, 250 /im of diameter), 
in order to eliminate the non-gaussian components and 
then it is collimated by a system of lenses to a diame- 
ter oi Wp = 1.25 mm. After the crystal, the pump is 
stopped by a couple of UV mirrors, transparent to the 
visible (T=98% at 710 nm), and by a low frequency-pass 
filter (T=95% at 710 nm). The down converted beams 



(signal and idler) are separated in polarization by two 
polarizers (T=97%) and finally the far field is imaged 
by a CCD camera [D We used a 1340X400 CCD array, 
Princeton Pixis:400BR (pixel size of 20 /im), with high 
quantum efficiency (around 80%) and low noise in the 
read out process (A = 4 electrons/pixel). The CCD ex- 
posure time is set by a mechanical shutter to 90 ms, thus 
each image cached in this time window corresponds to 
the PDC emission generated by a single shot of the laser. 
The far field is observed at the focal plane of the lens 
(/ = 10 cm) in a / — / optical configuration, ensuring 
that we image the Fourier transform of the crystal exit 
surface. Therefore, a single transverse mode q of the 
down converted photons is associated to a single point 
X = (A//27r)q of the image. 

We expect a strict correlation in the number of photons 
Ni and Ns detected at any couple of symmetric positions 
with respect to the center of symmetry, providing the de- 
tection areas are larger than the typical coherence area 
|23| . The coherence area Acoh^ i-e. the uncertainty on 
the relative position of the correlated photons in the far 
field, depends on the pump angular bandwidth [ij, and 
on the parametric gain [23, l28|. In our framework the 
detectors are the pixels of the CCD array and therefore 
the single pixel should be of the same order of magnitude 
of the coherence area or bigger. By contrary, if the pixel 
is smaller, the photons detected at the pixel in the po- 
sition X will have the correlated ones spread on several 
pixels around the position — x. This spoils the pixel-pixel 
correlation. 

For this reason, and in order to reduce the contribution 
of the read noise of the CCD, it is convenient to perform 
an hardware binning of the physical pixels. It consists in 
grouping the physical pixels in squared blocks NxN, each 
of them being processed by the CCD electronics as single 
" superpixel" (the photons collected by the superpixel are 
the sum of the photons of each pixel, whereas the read 
noise is just slightly increased with respect the one of the 
single pixel). As described in [23], following the previous 
argumentation, the size of the superpixel Apix must be 
at least of the order of the coherence area, although the 
ideal condition is Apix/Acoh ^ 1. In the present paper 
we consider three different superpixel sizes: 12x12, 24x24 
and 32x32, being the measured coherence area Acoh ^ 
(6pixel)2 = (120/im)2. 

The number of temporal modes collected in one shot 
image is Mtemp = {Tpump/Tcoh) ^ 5-10^, estimat- 
ing the coherence time Tcoh of the PDC process around 
one picosecond [V\. The number of spatial modes col- 
lected by the superpixel is M spatial ^ Apix/Acoh de- 
pends only on its size, since Acoh is fixed by the pump 
transverse size. From the total number of modes A4 = 
M^temp ^ M-spatiai ^ 10"^ ^ 10^ and the mean photon 
number per pixel one expect the excess noise of the multi- 
thermal statistics En = {Ns) /A4 '^0.1. It can be seen 
as the number of photons per mode, and here is fixed, 
since we keep fixed (aside unwanted fluctuation pulse- 
to-pulse) the power of the laser. Therefore the multi- 
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FIG. 2: Degree of Correlation and Fano factor versus the num- 
ber of photons over a sample of 400 frames analyzed. Fig2(a) 
presents the data for binning 12x12, i.e. a superpixel with size 
(240/im)^ and regions containing 320 superpixels. The mean 
value of correlation degree for this range of photon number 
is cr = 0.681 ± 0.003 and the standard deviation of the sam- 
ple is 0.06. The mean value of the Fano factor is F = 1.17 
and the standard deviation of the population is 0.13. Fig2(b) 
presents the same analysis for binning 24x24 (superpixel size 
(480 /j^m)^) and regions containing 56 superpixel. The mean 
value of correlation degree is a = 0.498 =b 0.004 and the stan- 
dard deviation 0.08. Fano factor has mean value F — 1.26 
and standard deviation of 0.35. 



thermal statistics of the source in the single channel is 
quasi Poissonian, and this is confirmed by the Fano fac- 
tors F = I ^ En = {5'^Ni) / {Ni) presented in Fig El 

The first step of the measurement consists in evalu- 
ating the degree of correlation between pixel pairs be- 
longing to large areas of PDC emission, without object. 
Chosen two correlated regions of PDC in the far field As 
and Ai of about {3mm)'^ size, we evaluate the degree of 
correlation in Eq. ^ by performing spatial averages over 
the ensemble of correlated pixels belonging to the two re- 
gions. The measurement procedure is described in |23| . 
Fig|2] shows the results obtained for a sample of 400 im- 
ages (each one corresponding to a pump shot and to a 



point in the graphs) for the two binning size. In both 
cases the degree of correlation is weh below the SNL. As 
expected, the result for binning 24x24 is better than for 
12x12, and the last is better then the one reported in [23|, 
anyway. Thus, we are in the regime a < 1 needed for ap- 
preciating an improvement of the SSNQI with respect to 
the classical differential scheme. 



IV. EXPERIMENTAL IMPERFECTION AND 
THEIR CORRECTION 

In this section we aim to illustrate the experimental 
challenges for achieving a strong local sub-shot-noise cor- 
relation with a large number of detected photons over a 
large number of pixels at the same time. Actually work- 
ing with high photon number is demanded for reducing 
the effect of the background. In fact the background 
(electronic noise of the CCD and straylight) can be cor- 
rected a posteriori on average, but the single image re- 
mains affected by this noise. Therefore we need to work 
in a condition of {S'^B)/{Ni -\- Ns) "C 1, where we indi- 
cate as B the number of background photons. The main 
sources of experimental excess noise are the following: 

1) the intensity gradients in the regions Ag/i of interest, 

2) the error in the determination of the center of sym- 
metry of the spatial correlation, 

3) the background noise due to the straylight and the 
read-out noise of the CCD. 

We analyze them in detail in the following subsections. 



A. Intensity gradients 

Up to now, we have considered an ideal situation in 
which the quantum efficiency rj is constant and equal in 
the two optical paths. 

From an experimental point of view, we had to cope 
with several sources of spatial dishomogeneity that con- 
tribute to the experimental excess noise and that can hide 
the quantum correlation, especially when large areas of 
PDC emission around degeneracy (710nm) are consid- 
ered. 

In the angular bandwidth corresponding to an area 
of (3 mm)^ close to collinear direction, lies a frequency 
bandwidth of about 10-20 nm. The response of the fil- 
ters and the CCD quantum efficiency have a frequency 
dependence even within this small range. This generates 
a spatial gradient in the intensity pattern. For instance 
the CCD sensitivity, that has a maximum around 800 
nm, varies of about 3% in our range of wavelengths. The 
bandpass filter (see Fig(T]) has a maximum of transmit- 
tance of 95% at 710 nm, decreasing of few parts per- 
cent at the border of our frequency region. These ef- 
fects are even more complicated by the non-linear de- 
pendence of the frequency on the emission angle in the 
phase-matching function. Even the optical transmission 
at different angles can be influenced by the polarization. 



Moreover, the spatial statistics of the detected light is 
affected by spatial imperfections of the optics, by the 
presence of dust and by the inhomogeneity of efficiency 
of the CCD pixels (a minor contribution). All these ef- 
fects contribute to an increasing of the excess noise and, 
consequently, a deterioration of the noise reduction fac- 
tor. 

We modeled it by introducing a spatial dependence of 
the quantum efficiency r]j^k^ with j=s,i and k identifying 
the symmetrical pixels in the position x/c and — x/e of the 
regions Ag and Ai respectively. The expectation values 
of the number of photons on the pixel k is: 



(Nj^k) =r]j,k{n) 



(4) 



with (n) the number of photon impinging a pixel in each 
frame. The associated variance and covariance are 129| 



(S^Nj^k) = VjMn)+vlk{nf/M 



(5) 



{5Ns,kSNi,k) = Vs,kViA{n) + {nf/M). (6) 

Using Eq. (g]), Eq. (0 and Eq. (P in Eq. (HJ we can 
calculate the theoretical value of the NRF: 
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where 7^+^^ = ^^'^2^''^ ^^ ^^^ mean value of the quantum 
efficiencies of signal and idler arm and r]-^k = Vs,k — Vi,k 
is their difference. It is worth to observe that the degree 
of correlation is deteriorated by the unbalancing in the 
quantum efficiencies r]-^k on the two paths. Due to the 
high pulse-to-pulse instability of the Q-switch laser, in 
our experiment it is convenient to estimate the NRF in 
the single frame obtained by a single laser shot. A single 
frame contains several pairwise spatial modes of the ra- 
diation, since the PDC is spatially broadband [16]; there- 
fore, pixels belonging to regions Ai and As are pairwise 
correlated in the number of photons, and the pairs can be 
considered as independent if the condition Apix ^ Acoh 
is fulfilled. We estimate NRF as defined in Eq.([T]), where 
the quantum mean values are evaluated by averages over 
the ensemble of the pixel pairs. The operator associated 
to the estimation of the mean photon number is 



1 ^ ^ 
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where 72. is the number of pixels in the region. According 
to Eq. dH, the expected value oi Af becomes 



1 ^ 
{Afj) = ^Y^njAn) = Vj{n)- 



(9) 
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Similarly, the operator for the variance of the photon 
number is defined as: 



^1 



lj2iN^^,-{J^^)f 



n 



(10) 



k=l 



Using Eq. (J5j), its expected value is {S'j) = r]j{n) + 

{rfj + y{Vj)/M) {n)'^/M + V{r]j){h)'^ (here we exclude 
correlation between the efficiency distribution in the two 
regions). We observe that performing the spatial statis- 
tics in order to estimate the theoretical quantities intro- 
duces a term of excess noise depending on the spatial 
variance V{r]j) of the detection efficiency in the region 
Aj. At the same time, we can define the operator asso- 
ciated to the estimation of covariances as 
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function of binning size and their values corrected by the 
flat field and for background noise. From the figures we 
notice that the effect of disuniformity, V{r]j)^ dominates 
and the flat field correction becomes fundamental. In 
fact, an increasing of the superpixel size corresponds to 
an increment of the number of photons and, consequently, 
of the excess noise generated by V{r]j) in Eq. (fT2|) . It is 
important to stress that this compensation is the same 
for all the images, since the transmission-detection ef- 
ficiency distribution is always the same. Therefore, this 
technique can be applied successfully even in the imaging 
scheme, when an object is inserted. 



k=l 



According to ([6]) and to (|9]),its expected value is (S^^) 

(n)_2 
M 

above, we can now estimate the experimental noise re- 
duction factor: 
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Comparing the expected value of the NRF a^^^ esti- 
mated in the single frame with the theoretical value ([7]) 
that refers to the single pixel pair, a new term of excess 
noise depending on the spatial variance of the efficiency 
appears. Although the variance is usually rather small, 
this term can be relevant (even the dominant one) when 
the number of photons per pixel becomes as high as we 
deal in our working condition, and needs to be addressed. 

This additional term can be corrected a posteriori by 
compensating the channel losses represented by rjj^k- Fo- 
cusing, as usual, on the two regions of interests, the com- 
pensation consists in multiplying each matrix-image by 
a sort of "flat field" matrix gj^k- This is obtained by the 

sum of Q single shot images Fj^k = X1q=i ^jk iioi"i^al- 
ized for the mean value in the region j as 
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The flat field correction corresponds, in the theoretical 
model, to the substitution Nj^j. -^ 9j,kJ^j,k in the expres- 
sions (|8|), (p!Q|) and (fTT|) . Consequently, applying the flat 
field compensation as in Eq. (fT3]) , the experimental NRF 
becomes: 



1 rjl (n) 



(14) 



that is a good estimation of the average NRF ([7|) over the 
two large regions. Thus, the application of this compen- 
sation for the losses disuniformity allows in principle to 
retrieve the expected quantum mean value of the NRF. 
Fig. [3] and Fig. 2] report the NRF and the Fano factor in 



B. Determination of the center of symmetry and 
sizing of the superpixel 

We observe that the theoretical prediction a = 1 — rj 
relies on the assumption that each pixel of As detects all 
the spatial modes correlated to the ones collected by the 
corresponding pixel in Ai^ and viceversa. Otherwise, the 
presence of uncorrelated modes in the two regions would 
not provide a complete cancelation of the excess noise by 
the subtraction [3l| . In particular it can be demonstrated 
that 
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with Mc/u the number of spatial corre- 
lated/uncorrellated modes detected by a pair of 
symmetric pixels. When the center of symmetry (CS) is 
estimated exactly, Mu = ^ and one retrieves the theoret- 
ical lower bound cf' = a. Therefore, experimental precise 
positioning of the regions and accurate determination 
of the CS is fundamental. In |32| it is demonstrated by 
realistic simulation that an high level of correlation is 
obtained by finding the CS with sub-pixel resolution. In 
the present work the CS used in a measurement with a 
certain NxN binning grid (N=12 or N=24), is evaluated 
previously by using a more resolved grid (typically 
2x2) as follows. Considering the region As we find the 
proper position of the correlated region Ai evaluating 
the noise reduction factor <j as a function of the position 
of Ai (moved by integer steps). The optimal position 
corresponds to a dip in the value of a. Finally, taking 
two correlated pixels Xs and Xi , the center of symmetry 
is Xc = i^i^^, A different approach is based on a 
procedure of NRF optimization by micro-positioning of 
the CCD array [ssj. 

However, for a given error in the determination of the 
CS, increasing Apix reduces the ratio between the uncor- 
related and the correlated modes. Therefore, an appro- 
priate determination of the minimum size of the super- 
pixel can reduce this bias under an accepted threshold. 
Fig. [3] shows clearly a strong dependence of the NRF 
(after flat field correction) from the pixel size. This be- 
haviour is well explained by a geometrical model based 
on Eq. ([15]). 



C. Background noise 

Another source of deterioration of the spatial correla- 
tions is the background noise. The two main sources of 
background noise are the electronic noise in the read-out 
and digitalization process and the straylight. The elec- 
tronic noise of the CCD is A = 4 photoelectron/pixel. 
This contribution can be strongly reduced in terms of 
signal to noise ratio when binning is performed. Stray- 
light light is mainly due to the broadband fluorescence 
of the crystal, mirrors and filters hit by the strong UV 
pulse of the laser and to the residual UV light itself. The 
background noise can be estimated independently by ac- 
quiring an image without PDC light, simply by rotating 
the pump polarization of 90 degrees. The experimental 
value of the NRF ^ can be corrected a posteriori taking 
into account the mean value of the background: 
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where (Bj) is the mean photon numbers of the area 
Aj of the background image and (S'^Bj) is its variance. 
The same correction can be applied to the Fano factors. 
From Fig. [3] we notice that the electronic noise is the 
main source of deterioration of sigma for low binning 
size (small Apix)-, while for larger binning the effect of 
disuniformity, Virjj)^ dominates. This can be explained 
considering that the read-out noise of the pixel (or su- 
perpixel) is a fixed value, independent on the binning. 
Nonetheless larger superpixels detect more photons, re- 
ducing the effect of read-out noise {{S'^Bj)/{Nj) <C 1). 
The same considerations are valid for the Fano factors 
(FigH]). The straylight collected per superpixel is ex- 
pected to have a Poissonian-like statistics and obviously 
increases with the binning together with the PDC light. 
Therefore, its contribution, and the relative correction in 
the graphs, should be independent on binning. 

We want to stress once again that, differently from 
the case of the fiat field correction, the background can 
be corrected on average, but the single image remains 
affected by its noise. Thus, one can not apply any cor- 
rection for the background noise to the SSNQI scheme. 



V. QUANTUM VS CLASSICAL IMAGING 

We insert a weak absorbing object in the signal beam 
region Rs. It is a thin (less than lOnm) titanium depo- 
sition of 9mm^ area on a glass, representing a tt. The 
SSNQI image of the object is obtained by subtracting 
from the direct image registered in the signal branch, the 
correlated intensity noise pattern registered in the idler 
branch, as it is depicted in the left side column of FigEfa) 
for binning 12x12 and Fig|5jb) for binning 24x24. In or- 
der to compare the results of the quantum scheme with 
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FIG. 3: The plot reports the value of the Noise Reduction Fac- 
tor as a function of binning size N. The red (continuous) line 
refers to raw data, without any noise correction. The corre- 
sponding data after the background subtraction are reported 
in the green (dashed) line. The red (dotted) line refers to data 
corrected with Flat Field and the black (dotted dashed) line 
data takes into account both Flat Field and background cor- 
rections. 



4.5 
4 

3.5 
3 



2.5 

2 

1.5 

1 
0.5 









/ 




Backg. Corn 






> 


^\ 






•J^ 1 f^^w 1 

1 




^.^-^ 


- 






-ield Corr. | 


^^=*s==^*'^ <;^ 


^'^^ — 


^^^ 






Flat Field + Backg Corr. 



10 15 

BINNING 



20 



25 



FIG. 4: The plot reports the value of the Fano factors as a 
function of the binning size N. The colors are according to 
Fig.E] 



the corresponding DCI scheme, we simulate the last one 
by subtracting from the direct image a noise pattern that 
is not correlated at the quantum level with the one con- 
taining the object (see right side of FigE] (a) and (b)). 
It is obtained simply by choosing a region belonging to 
the idler field that is shifted by at least one pixel with 
respect the quantum correlated one. In order to simulate 
a coherent illumination in both branches, the experimen- 
tal spatial noise in the region As and Ai should obeys to 
a Poissonian statistics, i.e. it should be shot noise lim- 



Direct 
Image 



Correlated 
Noise 



5 10 15 



i1600 
1400 





Correlation with the Reference noise free image 



s 




E 


10 
15 
20 


s 


1 



1600 
1 1500 
1 1400 



Reference 



Uncorrelated 
Noise 




• Direct Image 
SSNQI 

• DCi 








^S:':-: 



5 10 15 



Direct 4^J N^^°° 4 TT 





2 4 


6 8 




'^rm 




2 


■ 
















Noise 6 
8 






i 


T 









i9500 
9400 
9300 
9200 
9100 



2 4 6 8 

4 Ji JSn^ 
8H_»i ■ 



J 9300 
1 9200 
19100 



SP M 



FIG. 6: Correlation coefficient between a reference noise free 
image and the single shot frame over 400 frames. Fig. 4 
present the results for the 12x12 binning: the mean correla- 
tion coefficients are Cssqni = 0.71, Cdci = 0.60, C direct — 
0.75 respectively for the quantum, the differential classical, 
and the direct imaging. 
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FIG. 5: Experimental imaging using correlated properties of 
PDC quantum light, where the color-map is associated with 
the selected photon number as reported on the righthand side 
of each image. In (a) the pixel size is (240/im)^ and the 
absorption coefficient of the object is around a = 12%. In (b) 
we use less absorbing object with a = 5% and pixel size of 
(480/im)^ (see the text for discussion). 



ited. As we discussed above, in our working condition 
the low gain and the large number of collected modes 
limit the spatial excess noise of the single beams to small 
values (see Fano factor in Fig|2]and Figl4]). Furthermore, 
we report the images obtained by direct classical imag- 
ing, that in the present configuration represents the best 
possible classical scheme. By observing FigEl it is quite 
evident an improvement of the SSNQI image both with 
respect to the DCI and direct classical one. In the case 
of larger binning in particular, the object in the DCI im- 
age is even undistinguishable {a = 5%), whereas it is well 
perceivable in the SSNQI frame. There is always a trade- 
off between the need of enlarging the detection areas for 
more squeezing effect, and to preserve, at the same time 
the spatial resolution of the image. In fact this is just 



FIG. 7: Correlation coefficient between a reference noise free 
image and the single shot frame over 400 frames. Fig. 5 
present the results for the 24x24 binning: the mean cor- 
relation coefficients are Cssqni = 0.664, Cdci = 0.483, 
C direct = 0.641 respectively for the quantum, the differential 
classical, and the direct imaging, with statistical uncertainty 
smaller than 10~^ 



a technical problem, that in principle can be solved by 
using larger diameter of the pump laser and crystals. 

In order to quantify rigorously the improvement in the 
sensitivity of the SSNQI, we consider two different figures 
of merit. 



A. Correlation function 

The first one consists in creating a noise free reference 
image of the object, and to compare the SSNQI and the 



classical (DCI and direct) images against the reference 
one. Operatively, the noise free image can be prepared 
by averaging over a large number of direct images, so 
that almost all the noise, due to multi-thermal statis- 
tics of light as well as to the electronics of the digital- 
ization process, is washed away. The obtained reference 
images for the two binning conditions are reported in 
FiglSl The comparison of the SSNQI and the classical 
schemes against the reference is quantified by the square 
correlation coefficient 






(17) 



where j = SSNQI, DCI, XssQi{^) = N,{-^) - 7V;(x) 
and Xdci{^) = Ni{—^ + a) — A^^(x). a is the shift vec- 
tor and the prime (') indicates the number of photons 
detected in presence of the object. In principle the im- 
age that retains more detail of the object, inspite the 
deterioration due to the noise, has a larger correlation 
value Cj. For instance, for the images of FigEfb) we 
have CssNQi = 0.56, Cdci = 0.41, Cdirect = 0.49. 

The analysis has been performed on a sample of 400 
images and the results are reported in FigE] for binning 
12x12 and Fig |7] for binning 24x24. The correlation coef- 
ficient of the SSNQI (green marks) is always larger than 
the one for DCI (blu marks). In order to be completely 
fair one could take into account a sample of images that 
presents a null excess noise on average, i.e. a Fano fac- 
tor exactly equal to 1.0 on average. In fact we did, and 
the results do not differ from the behavior of the data 
presented in Fig [6] in the region of low photon number. 
As expected, for the larger superpixel size, the discrim- 
ination between the classical and the quantum imaging 
is more pronounced although more spread. This is re- 
lated with the degree of correlation that increases with 
the detection areas, as we have seen in Sec lIIII In the 
case of large pixel size, the SSNQI result is even better 
than the direct image obtained in the signal branch. This 
is in agreement with the value of the correlation degree 
that in the last case is smaller than 0.5 (see caption of 
Fig. [5fb)). This is the evidence that the method could 
becomes in the future a real tool for low noise imaging, 
whenever a small illumination level is needed. 



B. SNR 

The second figure of merit that we consider is strictly 
related to the theoretical Eq.(j3]). It consists in measur- 
ing the signal to noise ratio in the classical imaging and 
in the quantum imaging schemes and to compare them. 
For the following analysis it is convenient to insert in 
the signal beam an object with a uniform absorption 
a(x) = a ^ 5%, i.e. a square shaped deposition that 
can be considered having a constant thickness for our 
purposes. The prediction of Eq. (|3|) are tested on a 




a) 




b) 

FIG. 8: (a) R^^^ of the SNR in the SSNQI scheme and in the 
DCI scheme for binning 24x24. The blue marks represent the 
result that take in to account all the frames, without selection 
on the Fano factor, while the red marks are obtained taking 
into account only images that have Fano factors of F = 1±0.2 
(b) i?^^^^^* of the SNR in the SSNQI scheme and in the Direct 
scheme. The lines represent the theoretical predictions 



large number of frames (about 1200) for the case of bin- 
ning 24x24 and 32x32. We group the frames in classes 
with respect to their value of a, estimated for each im- 
age as described in Sec. IIIII Belonging to the j- class are 
the images with degree of correlation falling in a small 
band (bandwidth=0.1) with average (Jj and average Fano 
factor Fj . The estimation of the absorption spatial dis- 
tribution in the SSNQI scheme is obtained by following 
the Eq. ([2]). Experimentally, for the generic /c-th image 
belonging to the j-th class we perform the subtraction 

(7Vf^^(-x) - 7V^^^^(x))/(7Vf ^) = afc(x) between signal 

and idler correlated pattern and {N^ ^) = ^y^N\ (— x). 
Then the value of Qf(x) is obtained averaging over the 
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FIG. 9: (a) R^^^ of the SNR in the SSNQI scheme and in the 
DCI scheme for binning 32x32. The blue marks represent the 
result that take in to account all the frames, without selection 
on the Fano factor, while the red marks are obtained taking 
into account only images that have Fano factors of F = lzl=0.2 
(b) i?^^^^^* of the SNR in the SSNQI scheme and in the Direct 
scheme. The lines represent the theoretical predictions 



frames of the j-th class, i.e. 

^ AG 

a(x) = TT^^k{^) = Q^fe(x) 



(18) 



/c=l 



where J\fj is the number of images in the j-th class. 

The differential classical scheme with coherent beam 
can be simulated in our set-up by evaluating the absorp- 
tion as in the previous Eq. (p!8|) where the region of the 
idler is slightly shifted from the quantum correlated one 
Ni{—:si + a) (typically we take a = (1,0)). 

For each class, the SNR for the quantum and for 
the classical schem e s can b e estimated punctually as 
SNRj{x.) = a/c(x)/[J2Qf/j;(x)]"'^/^, where the overline in- 
dicates the average over the frames of the class j. 



At the same time we evaluate the factor R^^^{-x) = 

SNRf^^^/SNRf^^. Then, we perform also a spa- 
tial averages over x, to obtain a representative value of 
R?^^ for each class. The results are shown in Fig|8l^a) 
and Figini^a) against the correlation degree cfj for the 
binning 24x24 and 32x32 respectively: the SSNQI has, 
as expected, an advantage that increases with the level 
of correlation. Selecting only images with Fano factors 
very close to 1 for simulating a classical coherent illumi- 
nation (red data-set), the result are just slightly worse as 
one without Fano selection (blu data-set). However, for 
the best achieved values of the correlation degree in the 
quantum regime, we obtain an improvement of the SNR 
larger than than 70% compared with the classical imag- 
ing scheme. The results are in agreement with the theo- 
retical predictions obtained by the generalization of Eq. 
(|3]) with the substitutions E^"^^^ ^ Fj - 1 and a ^ aj. 

In Fig. [Htb) and FigEfb) we report the Ratio R of 
the SNR in the SSNQI scheme and in the direct classical 
scheme. When Fano values are selected to be ~ 1 (red 
data-set), allowing a proper simulation of a shot noise 
Hmited classical imaging, the results are in perfect agree- 
ment with what expected from theory. In particular for 
small a, we have R = ^==, i.e. the SSQI overcomes the 
direct classical imaging as soon as a < 0.5. 



VI. CONCLUSION 

We have presented the first experimental quantum 
imaging protocol aimed to improve the sensitivity of 
imaging techniques beyond the shot-noise-limit, that rep- 
resents the threshold between classical and quantum 
world. We reached a considerably high level of multi- 
mode spatial correlations using the properties of twin 
beams generated by the parametric-down-conversion pro- 
cess, and we used it for reducing the noise in the differ- 
ential imaging scheme. Our results represent the proof 
of principle of this new quantum technology. Anyway, 
some technical limitation toward practical application of 
the methods are still present in our setup. Firstly, total 
transmission of the optical channel should be improved 
in order to reach even higher quantum correlations. Fur- 
thermore, there is always a tradeoff between the need of 
enlarging the detection areas for more squeezing effect, 
and at the same time to preserve the spatial resolution of 
the image. In fact this is just a technical problem, that in 
principle can be solved by reducing the coherence area of 
the single spatial mode, in particular using larger diam- 
eter of the pump laser and crystals. A work for solving 
these problems in view of practical applications in now 
going on [26|. 
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